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Abstract 

We search for simultaneous baryon and lepton number violating decays of the meson. Specif- 
ically, we use 281 pb"^ of data taken on the ^(3770) resonance with CLEO-c detector at the CESR 
collider to look for decays — >■ pe+, pe^ , pe" and 5° — pe~ . We find no significant 

signals and set the following branching fraction upper limits: pe^{D^ pe'^) < 1.1 x 10~^ 

and L>° pe'iD^ ^ pe") < 1.0 x 10"^, both at 90% confidence level. 

PACS numbers: 13.20.Fc, 14.60.Cd 
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I. INTRODUCTION 



Various Grand Unified Theories (GUTs) [l| and many Standard Model (SM) extensions 
such as superstring models [2] and supersymmetry (SUSY) ^ predict baryon number vio- 
lation, and as a consequence nucleons can have finite, if long, lifetimes. However nucleon 
decay has not yet been observed [4]. In all these theories baryon [B) and lepton (L) number 
violations are allowed but the difference A{B — L) = is conserved. A higher generation 
SUSY model [5] predicts decay modes having such B and L violating decays for r leptons 
and for D and B mesons. The search for such r decays has been performed jl, 0], but decays 
of heavy quarks have not previously been investigated. 



In this paper we describe a search for the D meson decay channels D° — » pe^, 



pe~ 



— > pe~ and — >• pe~ . Such decays simultaneously violate B and L but conserve 
A(5 — L). Several models of proton decay, e.g. in GUT, superstrings and SUSY as described 
above can be augmented to provide predictions on possible decay mechanisms. 

In SU(5) theory, protons can decay into several modes; one of them is p — e^Ti^ . Biswal 
et al. [8] suggested five different decay diagrams. The decays are mediated by heavy hy- 
pothetical gauge bosons called X and Y . The X and Y bosons have electric charge |e 
and |e and couple a quark to a lepton, hence they are sometimes called "lepto-quarks." 
Figs. [H^a) and (c) show two of these possibilities that proceed via the s-channel. Fig. [It^b) 
is an analogous decay diagram for — > pe"*", where the mediator is a F boson. Here we 
take the coupling e^Yu as shown in Fig. [TJ^a) and introduce a coupling cYd replacing a u 
with a c in the t-channel version of Fig. [I](a). Similarly, Fig. [I](d) shows another analogous 
decay diagram for D*^ —>■ pe^ with an X boson as the mediator; here we take the coupling 
of e^Xd from Fig. [T]^c) and use the coupling cXu by replacing a u with a c in the t-channel 
version of Fig. [I](c). The spectator in both decay diagrams is u. No tree level diagrams allow 
— >• pe~ in SU(5). However, a decay model can be constructed using higher order dia- 
grams. Arnowitt and Nath also predict proton decay in an i?-parity violating [9;] superstring 





(a) 



u 



d d 
(c) 
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FIG. 1: (a) and (c) are s-channel decay diagrams of p 



described by SU(5) theory and (b) 



and (d) are decay diagrams of pe'^ based on analogous couplings. 



based model that can also accommodate pe+ decay [1 
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II. DATA SAMPLE, SIGNAL SELECTION AND RECONSTRUCTION METHOD 



(-) _ _ [— , 

We study the decays of pe~^ and pe using the CLEO-c detector [ll|]. We do not 

assume that the two modes, —>■ pe~^ and pe~ are charge-parity (CP) conserved 

(-) 

decays. When we refer to B{D^ pe^) we mean either B{D^ pe^) or B{D^ — > pe^). 

(-) 

Likewise, B{D^ pe~) is shorthand for either B{D^ P^ ) or B{D^ — * P^ )- 

The CLEO-c detector consists of a CsI(Tl) electromagnetic calorimeter, an inner vertex 
drift chamber, a central drift chamber, and a ring imaging Cherenkov (RICH) detector inside 
a superconducting solenoid magnet providing a 1.0 T magnetic field. In this study we use 281 
pb~^ of CLEO-c data produced in e+ e~ collisions and recorded at the ^(3770) resonance. 
At this energy, the events consist of a mixture of D^D^, D^D^ and e^e~ qq{q = u, d, s) 
continuum events with a small number of r^r~ and '~filj{2S) events. 

{-) 

We examine all the recorded events and look for candidates corresponding to — >■ 
pe^ and pe~ . The selection criteria for charged tracks are similar to that described in 
12l|. except that the momenta are required to be in the range from 50 MeV/c to 2 GeV/c. 



Moreover, we require that the polar angles, that the p{p) and e~{e~^) subtend with respect 
to the beam axis are required to satisfy | cos6\ < 0.9. Protons are identified using only the 
energy loss information [dE/dx) from the tracking chambers, since the kinematic limit of 
their momentum (900 MeV/c) is below threshold for radiation in the RICH detector. On the 
other hand, we do use the RICH, in combination with dE/dx, to aid in identification and 
elimination of kaons when the momentum is above 700 MeV/c, which is sufficiently above 



the RICH kaon radiation threshold. The specific requirements are discussed in Ref. [12 
Defining ap as the difference between the expected ionization loss for a proton and the 
measured loss divided by the measurement error, with analogous definitions for vr, K and e, 
we require \ap\ < 2.5, \(7t,\ > 3, \aK\ > 3 and — cXg < 0. We find that, for the momentum 
range of 0.5 to 0.9 GeV/c the proton identification efficiency is 98% and the probability that 
a pion (kaon) is misidentified as a proton is 0.9% (1.6%). 

Electrons (positrons) are selected as in Ref. [lH , with the additional criterion that we veto 
any candidate which passes the antiproton (proton) selection. The electron identification 

efficiency is 95%, with pion and kaon fake rates ~ 1%. 

(-) {-) 

We reconstruct candidates for both —>■ pe^ and pe~ modes separately. We 

evaluate the difference between the beam energy and the sum of the electron and proton 
energies (Ai?), and require \^E\ to be within two standard deviations (ctae = 5.3 MeV) of 
zero. For selected events, we compute the beam-constrained mass [3], defined as: 



where -Ebeam is the beam energy and pj represents the momenta of each final state particle. 
A signal would appear as a peak at the mass [3]. 



III. BACKGROUNDS AND SIGNAL SIMULATIONS 

Monte Carlo (MC) simulations are used to understand the response of the CLEO-c detec- 
tor, to characterize and estimate the possible backgrounds, and to determine efficiencies of 
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the reconstructed D° and decay modes. In each case e^e~ — >■ '?/'(3770) DD events are 
generated with the EVTGEN program fl5*|, and the response of the detector to the daugh- 
ters of the DD decays is simulated with GEANT 
simulation of initial state radiation (ISR) events, 
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I.e. 



. The EVTGEN program includes 
events in which the e+ and/or e~ 



17 is used to simulate 



radiates a photon before the annihilation. The program PHOTOS 
final state radiation (FSR). We use two types of MC events: 

• continuum MC events, in which e^e~ annihilations into uu, dd and ss quark pairs are 
simulated. It also includes the photon radiation by the initial state quarks. 

• signal MC events, in which either the D^ or the D° always decays in one of the 
two modes measured in this analysis while the other D^ or respectively, decays 
generically. 

The decay of D mesons into baryon pairs is kinematically forbidden, and so in the SM 
any real proton detected must be from a continuum event. Our largest source of potential 
background is the combination of a real proton from such an event with an electron from 
a photon pair conversion. We studied this background using a continuum MC simulation 
with five times the luminosity of our data sample. In Fig. [2] we plot the cos distribution, 
where is the angle between the e~ and any other candidate. All selection requirements 
are applied, except that we relax the /S.E requirement to ±4(T, and accept candidates in the 
broader Mbc range between 1.83 and 1.89 GeV. A clear excess near cos0 = 1 is observed. 
We remove these events by requiring cos0 < 0.73, which removes 71% of the background 

with a 3.4% loss in signal efficiency. 

(-) 

We determine the D^ signal line shape parameters and detection efficiencies using a signal 
MC sample for each mode. The Mbc distributions are shown in Fig. [31 We describe the 
signal shape using a Crystal Ball function [l8|, which has the form: 



/(MbclMcCTAf, 



,a,n) 



A exp 



Mb 



A 



here A 



'be 

-1 . 



for Mbc < Md - olctm^, 
for Mbc > Md - aau,,. 



n 1 ^ 
a n—1 



+ 



1 



where A is an overall normalization, Md is the D^ mass i^, (Tm^^ is the mass resolution, 
and n and a are parameters governing the shape of the high mass tail. This high mass tail 
results from initial state radiation from the e~ and/or beams. In each fit, the parameters 
are determined by a binned maximum likelihood fit and their values are fixed in fits to data, 
with the exception of n. The fits are highly insensitive to the precise value of n, which is 
fixed to 7.0 throughout the analysis. 

(-) (-) 

From the reconstructed yields, we determine signal efficiencies for D^ pe^ and D^ — >■ 
pe- of (59.1±0.5)% and (59.4±0.5)%, respectively. 



IV. RESULTS FROM DATA 

The Mbc distribution of events passing all selection criteria in data are shown in Fig. H] 

(-) (-) 

for D^ pe^ and D^ — > pe~ modes separately. The background shape is parameterized by 
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FIG. 2: Distributions of cos0, the angle between e+e candidates, as discussed in the text, for (a) 



(-) 



continuum MC, and (b) data. The dotted histograms show cases where the e"*" is from a — > pe"' 

pe~ 



{-) 





candidate, and the sohd histograms correspond to cases where the e comes from a D 
candidate. Events to the right of the arrows are ehminated. 



an ARGUS threshold function [19], which has the form 



1 - 



be 



beam 



exp 



5 ^1 - 
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(2) 



Here, K is an overall normalization, and the other parameters, Ebeam and S, govern the shape 

of the distribution; Eheam is the beam energy and S" is a scale factor for the exponential. We 

(-) (-) 

fit the Mbc distributions of the individual modes — * pe~^ and pe^ with fixed signal 

shape parameters (from the signal MC) and fix Eheam = 1-8865 GeV, but fioat parameters 
K and S in the background function. The fits are shown in Fig. 111(a) and Fig. Hl^b). In both 
cases, the fit yield is zero and upper limits will be computed. 



V. SYSTEMATIC ERRORS 

We consider a number of systematic errors. We assign ±0.7% systematic error for finding 
each charged track, hence ±1.4% for both tracks. For electron identification we assign 



±1% error [13|. The proton identification uncertainty is ±1%, and was evaluated at higher 
beam energies on/near the Upsilon resonances by comparing the efficiency for identifying the 
proton in A — > pn~ decays in data and Monte Carlo simulation. Additional cross-checks were 
performed at center of mass energy, Eqm = 3770 MeV that showed consistent performance of 
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1.830 1.840 1.850 1.860 1.870 1.880 1.890 



M,, (GeV) 

(-) (-) 
FIG. 3: Mbc distributions for (a) — > pe~^ , and (b) for — > pe~ from signal MC, fitted with 

Crystal Ball functions. 



the particle identification over these running periods. Thus the overall particle identification 
uncertainty is ±2%. 

To estimate the systematic error arising from the AE cut, we compare signal yields using 
the nominal AE cut and a wide AE cut of ±100 MeV for the kinematically similar — »• 
K^Ti^ decay. The fractional decrease between the nominal and wide AE cuts are 9.02 ± 
0.34% and 8.93 ± 0.14% for data and MC simulation, respectively. The difference is 0.09 ± 
0.37% and therefore we assign a systematic uncertainty of ±0.4% to account for possible 
mismodelling of this quantity. The selection of cos0 < 0.73 reduces the efficiency by only 
3.4%, and we assign a ±1% uncertainty to the efficiency due to this cut. The uncertainties 
in the background shape due to the threshold Ebeam are determined by calculating the 
differences in the 90% confidence level, (C.L) upper limit yields between the nominal fit 
and a fit with Eheam shifted by ±0.5 MeV. Then the differences in the upper limits at the 
90% confidence level (C.L.) were taken to be the systematic errors from this source, which 
we estimate as ± 1%. Uncertainties due to the signal shape parameters were found to be 
neghgible. We also sought possible uncertainties due to differences in the veto efficiencies 
between data and simulation for kaons faking protons and antiprotons faking electrons and 
similarly for the charge conjugates. The differences were negligible. Finite MC statistics 
also introduces a 0.8% systematic error. The systematic errors are summarized in Table HI 

VI. UPPER LIMITS OF BRANCHING FRACTION 

The likelihood distributions as a function of the assumed yields are shown in Fig. [5] for 
(a) the sum of possible —>■ pe~ and —>■ pe~ yields, and for (b) the sum of possible 
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1.830 1.840 1.850 1.860 1.870 1.880 1.890 



M,, (GeV) 

(-) (-) 
FIG. 4: Mbc distributions for (a) pe~ and (b) for pe~^ from data shown by solid 

histograms. The curves are the fits as described in text. 



TABLE I: Systematic Uncertainties 

Sources of errors error(ib%) 

Tracking 1.4 

Particle identification 2 

cut 0.4 

COS0 cut 1 

Background shape 1 
Relative statistical error from signal MC 0.8 

Total in quadrature 3.0 



— >• pe~^ and — » pe~^ yields. We determine the upper limits of branching ratios by 
integrating the likelihood function to include 90% of the probability. We find 90% confidence 
level (C.L.) upper limits of 6.40 and 6.00 events, respectively. We compute the upper hmits 
on the branching fractions using: 



Here, iV^oso = (1.031 ± 0.008 ± 0.013) x 10^ is the number D'^D^ events at the ^(3770), 
where the first error is statistical and the second is due to systematics [20], is the 90% 
C.L. upper limit and e_ is the signal MC efficiency, reduced by one standard deviation. We 
determine an upper limit for the sum B{D^ — > pe~^) + B{D^ —>■ pe~^). We interpret this as a 
conservative upper limit on B{D^ pe~^) or B{D^ pe^). A similar interpretation is used 
(-) 

for B{D^ —>■ pe~). The calculated upper limits with and without the systematic errors are 
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(-) 



shown in Table [TTl In particular, we find B{D^ 
1.0 X 10-^ both at 90% C.L. 



(-) 



pe+) < 1.1 X 10-5 and ^ pe") < 




(-) (-) 
FIG. 5: Fit Likelihood plots versus the yield N for (a) for — > pe"*" and (b) for pe~ from 

data. In each plot, the vertical line shows the value of below which 90% of the total area lies. 



TABLE II: Results from fits to the My,c distributions and the resulting upper limits on branching 
fractions for both of the modes. 







(-) 


(-) 






L>o ^ pe+ 


L>o pe" 


Upper limit 


on N 


6.42 


5.94 


Upper limit on N(includin| 


5 systematic errors) 


6.61 


6.12 


Upper limit 


on B 


< 1.1 X 10-5 


< 1.0 X 10-5 



VII. CONCLUSIONS 

(-) (-) 

We have searched for the B and L violating decays pe^ and pe~ and 

find no evidence of these decays. We obtain branching fraction upper limits of B{D^ — >■ 
pe+)[i3(L>° ^ pe+)] < 1.1 X 10-5 and ^ pe-){B{D'^ pe-)] < 1.0 x 10-5, both at 

90% C.L. Using these limits, and the Ufetime, r^o = (410.1 ± 1.5) fs we compute the 
partial widths (Fj = Bi/tdo) to be: 

(-) (-) 
r(L'° pe+) < 2.8 X 10^ s-^ and F(D° pe') < 2.5 x 10^ s'K (4) 

These decay width limits provide less stringent constraints on new physics interactions than, 
for instance, proton decay experiments. However, no previous searches have investigated the 
possibility of charmed mesons violating B and L. These limits do not violate the predictions 
of higher generation models, which predicts pl^) ~ 10-^^ [5|. 
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